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Abstract 

The problem of black hole production in transplanckian particle collisions is revisited, in the 
context of large extra dimensions scenarios of TeV-scale gravity. The validity of the standard 
description of this process (two colliding Aichelburg-Sexl shock waves in classical Einstein gravity) 
is questioned. It is observed that the classical spacetime has large curvature along the transverse 
collision plane, as signaled by the curvature invariant {R^uXa)'^- Thus quantum gravity effects, and 
in particular higher curvature corrections to the Einstein gravity, cannot be ignored. To give a 
specific example of what may happen, the collision is re-analyzed in the Einstein-Lanczos-Lovelock 
gravity theory, which modifies the Einstein-Hilbert Lagrangian by adding a particular 'Gauss- 
Bonnet' combination of curvature squared terms. The analysis uses a series of approximations, 
which reduce the field equations to a tractable second order nonlinear PDE of the Monge- Ampere 
type. It is found that the resulting spacetime is significantly different from the pure Einstein case 
in the future of the transverse collision plane. These considerations cast serious doubts on the 
geometric cross section estimate, which is based on the classical Einstein gravity description of the 
black hole production process. 

PACS numbers: 04.70.-s, 04.50.+h, ll.lOKk 



I. INTRODUCTION 



All present-day macroscopic experimental data are consistent with gravity being described 
by the 4-dimensional Einstein-Hilbert action 

S = ^ / Ry/^d^X (1) 



In the microworld no gravitational effects have been observed, and indeed the 4-dimensional 
(4d) Newton constant G in (P) is so small that such effects would be totally negligible up to 
unreachable collider energies of M^f ~ 10^9 GeV (the 4d Planck scale). 

The idea of large extra dimensions allows to imagine a world in which gravity is strong 
enough to play a role in elementary particle collisions at accessible energies. In this world 
Einstein gravity propagates in D > A dimensions, out of which D — A are curled up in a 
compact manifold of size i ^ l/Mp^ . At distances ^ i such gravity would be described 
by an effective 4d action of the form (0). To get the effective 4d Newton constant right, we 
have to set the D-dimensional Planck scale at 

il4f)~M(i)(Mg)£)-^«M(i). (2) 

This means that such gravity will be much stronger than we are used to believe at short 
distances <^ i. 

Motivated by the hierarchy problem, Arkani-Hamed, Dimopoulos and Dvali proposed 
to push this idea to the extreme and lower Mpi all the way down to the electroweak scale 
Mew ~ 1 TeV. This proposal requires extra dimensions ranging from a mm to a fermi for 
6 < D < 11 {D = 5 leads to £ ~ 10^'^ cm and is excluded). To make this consistent with the 
4d Standard Model, we have to assume that all the other fields but gravity do not feel the 
extra dimensions. They must be confined to a 4d submanifold (brane) of the D-dimensional 
world (bulk). 

The simplest collider signature of such braneworld scenarios would be apparent en- 
ergy non-conservation due to produced gravitons escaping into the bulk. Recent Tevatron 
sea..cKe. of such ,„,.i„g e„e.,y event, fou.d ,>o .a*.caU. s,g,.ficaM effect and set a 
lower bound of ~ 0.6 TeV on the D-dimensional Planck scale (for a general review of signals 
and constraints, see ^]). Of course the LHC with its cm. energy of 14 TeV will be a much 
better probe of such phenomena. 
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In this paper I would like to contribute to the unfinished discussion of another possible 
.ignatwe of TeV.cale g.av.ty - black hole production in t.an^plancHa. colU.on. (see 
reviews j^). 

The common current opinion is that this process may be adequately described using 
classical general relativity, and that a single large black hole will form for a range of impact 
parameters (in the Z)-dimensional Planck units) 

b < RscUE) E'/^''~'\ (3) 

Here ^ 1 is the energy of colliding particles, and Rschw{E) is the Schwarzschild radius 
of the D-dimensional static black hole of mass E. If production cross section 

a^7i[Rs,UE)Y (4) 

based on this estimate (known as the geometric cross section) were true, the LHC would 
produce black holes at a rate ~1 Hz for Mpi = 1 TeV, becoming a black hole factory j^. 

In older studies of transplanckian collisions O] the authors usually carefully stated that 
in the range (jH)) a black hole is expected to form. Indeed, the eikonal expansion, normally 
used to analyze quantum gravity in the transplanckian regime, breaks down exactly in that 
range. It seems that the geometric cross section proposal was first put forward unequivocally 
by Banks and Fischler |7|. 

Later on, an appealing theoretical argument backing 0, (jH) was proposed by Eardley 
and Giddings 0]. In the first step of the argument the problem of black hole formation is 
analyzed in classical Einstein gravity using the closed trapped surface (GTS) method. A well- 
defined mathematical problem of finding a GTS in the spacetime formed by superposing two 
Aichelburg-Sexl shock waves is formulated and solved. For the range of impact parameters 
when a GTS is found, black hole formation is concluded by invoking the Gosmic Gensorship 
Gonjecture. 

In the second step the aim is to argue that quantum corrections to general relativity are 
not likely to modify the conclusions, because the classical spacetime has small curvature in 
the regions relevant for the trapped surface evolution. However, it is this second step where 
I disagree with Eardley and Giddings.^ 



^ It should be mentioned that the vaHdity of the geometric cross section was also challenged by Voloshin 
01 . Some of his criticisms were addressed in liof , others still remain unanswered. There does not seem 
to be an obvious connection between Voloshin's and my reasons for critique. 
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As we will see, curvature becomes large on the transverse plane at the moment when 
the particles pass each other, and in the future of this plane corrections to Einstein gravity 
cannot be ignored. In particular, if higher order curvature terms are present in the effective 
gravitational Lagrangian (and we have no reason to believe that they don't), they will 
become important precisely at this moment. 

To see a specific example of how this might happen, we will add to the Einstein-Hilbert 
Lagrangian a particular combination of curvature squared terms {R^uXaY — '^iRfiu)'^ + 
first considered in D dimensions by Lovelock [11], with a coefficient of natural magnitude 
~ 1 in Planck units. The technical reason why I choose this combination is that it leads to 
second order equations of motion. My primary goal in this paper is to furnish an example 
of how things may go wrong, and it is unlikely that any other combination would be better 
in this respect. 

In this Einstein-Lovelock theory we will then consider the simplest — zero impact pa- 
rameter — collision case. Existence of a CTS for such head-on collisions was shown long 
ago bv Pe,n.„.e Q. We will see that the curvatu. .,ua.ed te™. .deed beco.e ..po.ant 
and significantly modify the spacetime to the future of the transverse collision plane. 




FIG. 1: The longitudinal slice of the collision spacetime in the pure Einstein and the Einstein- 
Lovelock theories 

The result of this analysis is summarized in Fig. ^ In the Einstein case the post-collision 
spacetime is weakly curved (at large transverse radii). The metric in the future wedge F 
can be expanded in a power series in the light-cone coordinates m, f , starting with a linear 
term (xu + v. The Riemann tensor has 5-function singularities on the null planes m = and 
f = 0; the Ricci tensor is of course everywhere zero. That curvature becomes large can be 
seen from the invariant {R^yx^^o:. 5{u)5{v). 



4 



On the contrary, in the Einstein-Lovelock case the future wedge is spht into sectors 
n > 2, separated by planes on which the Riemann tensor has additional S- 
function singularities. Moreover, the Ricci tensor is also singular on these planes, and in 
general nonzero in between. In particular, the post-collision spacetime does not satisfy 
Einstein's equations. An additional feature is that the number of sectors n is arbitrary; the 
solution is not unique. 

The paper thus contains two main results — a general observation (curvature becomes 
large) and a specific example (analysis in the Einstein-Lovelock theory) — which show 
that classical Einstein gravity may not be applicable to study black hole production in 
transplanckian collisions. If so, the geometric cross section estimate (0)) becomes much less 
certain, having lost its most serious supporting argument. 

The plan of the paper is as follows. In Sec. II we will review the GTS argument. 
Readers familiar with it may go directly to Sec. Ill, where we will discuss the danger posed 
by regions of high curvature and see that curvature becomes large along the transverse plane 
at the collision moment. This conclusion is reached after the finite width of the shocks due 
to classical and quantum smearing is taken into account. 

The remaining part of the paper deals with the Einstein-Lovelock theory. In Sec. IV we 
will review the Lovelock Lagrangian and its higher order generalizations. I will also comment 
on the history of this Lagrangian and its alleged connection to string theory. 

In Sec. V we will write down a metric ansatz for the head-on collision, in a plane-wave 
approximation. Then we will derive equations of motion, in a small-gradient approximation, 
in the region of longitudinal light-cone coordinates |f | ^ 1. We will also see that the 
cubic and higher order Lovelock Lagrangians do not contribute in this approximation. 

In Sec. VI we will solve the resulting hyperbolic Monge-Ampere equation. First we will 
see that no usual-sense (smooth) solution exist, even when we consider the smeared out 
shocks as initial conditions. Then, in the limit of infinitely thin shocks we will find an 
infinite class of generalized (weak) solutions. We will discuss what this non- uniqueness may 
mean, and how it can possibly be resolved. 

In Sec. VII we will check that the small-gradient approximation made in Sec. V was 
justified. We will then use the found solutions of the Monge-Ampere equation to write down 
approximate spacetime metrics. I will explain the differences between the Einstein and 
Einstein-Lovelock post-collision spacetimes, and their negative implications for the GTS 
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argument. I will conclude in Sec. VIII with a summary of the obtained results. 

Notation. In the remainder of the paper we work in flat spacetime of D > 5 dimensions 
(although we quote some D = 4 results for comparison). This is legitimate, since the relevant 
gravitational dynamics happens in a region of ~ TeV^^ size, which is much smaller than 
the size of extra dimensions. We use Planck units, setting 8ttG = 1 in the D-dimensional 
Einstein-Hilbert action. Greek indices /x, z/, . . . run through all D coordinates, Latin i,j, . . . 
only through D — 2 transverse directions. Spacetime signature is (+,—,—,...), and the 
curvature tensor definitions are R'^^ap ~ ^a^^p ~ ■ ■ ■ and R^^ = R^'^j^^u ('Landau-Lifshitz 
timelike conventions'). 



II. CLOSED TRAPPED SURFACE ARGUMENT 

The basic ingredient in existing discussions of black hole formation in particle collisions 
is the gravitational field of a fast point particle moving along a straight line. In the limit of 
infinite boost '-f = E/m —>■ oo and fixed energy E, this field takes the form of a shock wave 

ds^ = dudv — ^(x) 5{u) du^ — dx^. (5) 

Here u = i— z^v = t + z, the particle is moving in the positive z direction, and x denotes 
the D — 2 transverse coordinates. 

Einstein's equations with the lightlike source 

nu = E6{u)6''~\x) (6) 

give the following equation for the shock wave profile $(x): 

-^V2$ = E5^-2(a;). (7) 

Thus we get^ (f = '^d-z is the volume of the unit D — 3 sphere) 

E 

$ = Inf [p = 4), (8) 

IT 



^ This result can also be derived by boosting the static Z?-dimensional Schwarzschild solution and simulta- 
neously scaling down the rest mass to keep the energy constant. This was the method originally used by 
Aichelburg and Sexl in D = 4 |13j. 
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Form (0) of the metric is unsuitable for analyzing the behavior of geodesies crossing the 
shock at M = 0, which is necessary for understanding the causal structure. For this we 
must switch to coordinates in which the metric is continuous. This is accomplished by the 



discontinuous coordinate transformation 
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u 

V 



x' + -Vi^{x)e{u) 



(10) 

(11) 

(12) 



(where 9 is the Heaviside step function). In the new coordinates the metric becomes 



du dv — HikHjk dx^dx^ 



Hij = S^j + -ViVj<!>{x)ue{u), 



(13) 
(14) 



and both geodesies and their tangents are continuous across the shock. Introducing polar 
coordinates in the transverse plane, this metric can be written as (see for D = 4) 



ds = dudv — 



1 + 



{D - 3)E 



1 2 



—u6{u) 



D-2 



dr' 



1 - 



E 



D-2 



ue(u) r'dQ\ 



(15) 



To discuss a collision, we add a second fast particle of the same energy moving opposite 
to the first one in the negative z direction at a transverse distance (impact parameter) b. 
Assuming that the particles pass each other at u = v = 0, the combined gravitational field 
outside the wedge u,v > can be obtained by superposing the metric (fT3|) with its mirror 
image: 



ds'^ = dudv — {HikHjk + HikHjk — 5ij) dx''dx\ 
H,, = 6,, + ^V^Vj^{x-{b,0,...,0))veiv). 



(16) 
(17) 



Such superposition is legal, because the excised region m, t> > (wedge F in Fig. (T]) is 
precisely the future light cone of the collision plane of the shocks. By causahty, outside this 
region the shocks will not be able to influence each other. 

The metric for u,v > must be found by solving the characteristic initial value problem 
for Einstein's equations. This task is complicated by the fact that the shocks after passing 
each other will focus and develop regions of high curvature. Shock propagation through 
these regions depends on their detailed structure jl^ . 



Nevertheless, we are interested whether this comphcated dynamics will lead to black hole 
formation. Since the complete metric is unknown, we must resort to indirect arguments. 
The idea j3[y] is to look for a closed trapped surface (CTS) in the part (fTH|) of the spacetime 
that we do know. CTS is defined as a closed {D — 2)-surface whose area decreases locally 



when propagated along the outer nul 



surface have positive convergence 16 1. 



normals. Equivalently, the outer null normals of such 



Existence of a CTS in a spacetime solving vacuum Einstein's equations implies presence 
of a singularity in the future. Assuming Cosmic Censorship, this singularity must be hidden 
behind a horizon, and we may conclude that a black hole will form. Moreover, the black 
hole horizon must lie outside the CTS. Using this information, one can get estimates of the 
horizon area and, via the Area Theorem, of the mass of the formed black hole. 




FIG. 2: CTS (|19|) . The trajectories of the colliding particles and the null planes containing the 
CTS are also shown. All the transverse directions but are suppressed. 

For the head-on collision (6 = 0) a CTS is easy to find jl^. Isj]. It lies in the union of the 
shock planes m = and f = and, when transformed to the m, x coordinates, consists of 
two fiat {D — 2) -dimensional disks of radii 



,l/(D-3) 



(18) 



centered at the collision point. In the u, v, x coordinates the CTS is glued out of two halves 
described by 

u = 0, V = $(ro) — r < rn; 

" (19) 

f = 0, u = $(ro) — ^{r), r < tq. 
It is symmetric under rotations of the transverse directions and the reflection z —z (see 
Fig. 12)). Strictly speaking, this surface is a marginal CTS, which means that the outer null 
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normals have zero convergence. (Such a surface is also called apparent horizon.) However, 
moving a small distance inside one can find a true CTS with negative convergence. 

For nonzero impact parameters {h > 0) the search of trapped surfaces was carried out 



analytically for D = 
Yoshino and Nambu 
impact parameters 



: by Eardley and Giddings |S] and numerically for 5 < D < 11 by 
17[. In every D existence of a marginal CTS was established for 



/ 2E \ 1/(^-3) 
6<..flsewW=c,(^^^-^^) , (20) 

where cd ~ 1-5 is a numerical constant weakly depending on D. This result of classical 
general relativity if the most serious argument in favor of the geometric cross section for the 
black hole production rate. 

III. HIGH CURVATURE AT u = v = 

The CTS argument of the previous section was valid for the Einstein gravity with its 
vacuum equation of motion R^^, = 0. This equation is used crucially in the Raychaudhuri 
equation describing how the convergence 6 evolves along the congruence of null normals i^: 

^ = R,J>^r + + 2(shear)2. (21) 

For Ricci-fiat spacetimes the first term vanishes. One then concludes that once 9 goes posi- 
tive, it grows monotonically and becomes infinite at finite affine distance A. This argument 
is the basis of the proof of singularity theorems jl6 |. 

Now suppose that long before the singularity forms, the evolving surface passes through 
a region of high curvature. If higher order curvature terms are present in the effective 
gravitational Lagrangian, they will become important and modify the field equations at this 
point. In general, R^y will no longer vanish, and one cannot make any conclusions about 
the dynamics of 9 across the region of high curvature. 

From the effective field theory point of view, it seems likely that all interactions al- 
lowed by symmetry should become significant in a theory at its natural scale (which is 
the D-dimensional Planck scale in our case). This is the main reason to believe that the 
gravitational Lagrangian will contain higher curvature terms, with coefficients of natural 
magnitude ~ 1. There are also other hints pointing in the same direction. For instance. 



since the Einstein gravity is non-renormalizable (see, e.g., review [1q\), higher curvature 
counterterms should be included to cancel infinities in loop diagrams. In string theory, 
which is a finite theory of quantum gravity a particular sequence of higher curvature 
terms with calculable coefficients appears in the effective action (the a'-expansion).^ 

It should be noted that singularity theorems of general relativity can also be formulated 
and proved in presence of matter sources satisfying various positivity assumptions, such as 
the null energy condition Tfj_iyi'^i'^ > 0. While such conditions are natural for classical matter, 



2l|. 



they can be violated by the renormalized energy-momentum tensor of quantum fields 
So it is unlikely that any useful generalization of singularity theorems exists when quantum 
effects are taken into account. 

These considerations make it clear that appearance of high curvature regions is problem- 
atic for the whole CTS argument. Eardley and Giddings tried to address this problem. In 
particular they noticed that the CTS they constructed lies in the planes of shocks. They 
proposed to bring it out of these planes by propagating it a small affine distance to the future 
along the outer null geodesies. This way, they said, one can get a CTS lying everywhere in 
the region of small curvature. The conclusion was that quantum gravity effects are unlikely 
to modify the result: the black hole will still form. 

I believe that this argument as it stands is incomplete: it misses the fact that curvature 
may become large on the part of the trapped surface where it crosses the transverse collision 
plane. 

First of all, the only nonzero components of the Riemann tensor of spacetime © are 

Rmuj = ^^i^Mx)5{u). (22) 

All curvature invariants (contractions of this tensor with itself and the metric) vanish iden- 
tically. Thus a shock plane by itself is not a region of high curvature.^ 

However, a problematic region does appear when we add a second shock. It is located at 



^ For completeness, we have to mention that Loop Quantum Gravity attempts to quantize general rela- 
tivity non-perturbatively, starting from the pure Einstein-Hilbert action. Since the program is unfinished, 
and in particular the classical limit is poorly understood, it is hard to say what form the gravitational 
field dynamics will eventually take. It seems likely that the effective action, as long as this concept is 
applicable, will have to contain higher curvature terms even in this theory. 

^ This also follows from the fact that the Aichelburg-Sexl wave can be obtained by boosting from manifestly 
low-curvature regions far away from a static _D-dimensional black hole. 
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u = V = 0, where the shocks colhde. At this moment we can form a nonvanishing curvature 
invariant of the coUision spacetime ()16|) : 

Rf.uXaR^"^'' OC RmujRvivj Cx6{u)6{v). (23) 

This equation seems to suggest that curvature becomes large (and even infinite) all along 
the transverse plane u = v = 0. However, it would be incorrect to jump to this conclusion 
too soon. The reason is that infinitely thin shocks are an idealization; in reality the shocks 
will have a finite width w. 

There are essentially two reasons for w > 0. The first, purely classical, reason is that in 
practice j = E/m is large but not infinite. Because of this, the shocks have width 

Wclass^r/j, (24) 

depending on the transverse distance r. This width comes out naturally when deriving the 



shock wave by boosting the static black hole spacetime |22l \14\. 

The second reason has quantum nature. In relativistic quantum theory coordinates of 
a particle at rest cannot be measured more precisely than its Compton wavelength ~ 1/m 
(see, e.g., js^) . For an ultrarelativistic particle, this limit becomes ~ 1/E. Thus the 
point particle picture used in the CTS argument should not be taken literally. In reality the 
particles should be thought to have finite size 1/E, which translates to a contribution to 
the shock width 

?i;q,an ~ 1/E. (25) 

We are interested in curvature at transverse radii r ~ ro relevant for the trapped surface 
formation and evolution. In practice, for such r we will always have Wdass <^ u^quan, because 
the masses of colliding particles are small compared to 1 TeV, and thus 7 is very large. (E.g. 
7 ~ 10^ at the LHC, using current-quark masses mu,d ~ 10 MeV.) Thus we can estimate 

W~Wquan~l/^- (26) 

Because of nonzero width, 5-functions in ()22|) . ()23p will be smeared out over an interval 
of length Au ~ w. Since the integral has to remain unity, the maximal attained value will 
be ~ l/ty ~ i?. Using this value in (j2Sl), we find 

max {R^^xa? = -2){D- 3) r^^g (27) 

r=ro 
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This formula involves a positive power of E, and a large prefactor (increasing from ~ 300 to 
~ 670 for L) = 5...11). 

We see that even after finite shock width is taken into account, curvature still becomes 
large in the region of transverse collision plane relevant for the black hole formation. This 
means that corrections to the Einstein gravity will become important at this moment. Their 
subsequent effect on dynamics is hard to predict in general. It is not a priori excluded that 
this effect will be transient, localized in a small shock-interaction region to the future of 
u = V = 0, after which the gravitational field will return to its Einsteinian value. In this 
case we could push the CTS through this small problematic region, similarly to Eardley and 
Giddings's proposal, and get to small curvature values where the Einstein gravity would 
again be applicable. (In the original proposal 8, p. 6] the part of the trapped surface near 
u = V = was to be left fixed.) However, I think that such localization is unlikely to 
happen. It is much more probable that the corrections will be significant up to the values 
u,v ~ 1, and even further. This would mean that the Einstein gravity alone essentially 
loses its predictive power in the future of the collision plane (a possibility also mentioned by 
Kancheh |2J]). 

In the remaining part of the paper I will give an example of how this scenario may be 
realized, using a particular modification of the Einstein gravity by curvature squared terms. 

IV. EINSTEIN-LOVELOCK GRAVITY 

We consider a modification of the D-dimensional Einstein gravity by curvature squared 
terms described by the Lagrangian 

C = + (28) 

jO^^^ = {Rf.uXaR^''^''-^Rf.uR^'' + R^)V^- (29) 

The new coupling k, is assumed to be ~ 1 on the grounds of naturalness. Lagrangian ()29|l 



has an interesting property first discovered by Lovelock it is the only combination of 
curvature squared terms leading to second order equations of motion for the metric. (A 
generic combination would produce fourth order equations.) 

For this reason theory (j28p has been often used in attempts to understand how higher 
curvature corrections modify the behavior of pure Einstein gravity. Thorough reviews of 
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the existing hterature can be found in j25| (apphcations to black holes) and [26[ (Kaluza- 
Klein scenarios; braneworlds) . Recently, theory was even discussed in connection with 
black holes produced in particle collisions: the authors of |27] argued that by studying their 
Hawking evaporation spectra one can measure k. However, in this paper we are interested 
how the Lovelock term manifests itself before rather than after the black hole formation, and 
whether it may in fact preclude this formation. 

Actually, Lovelock [ll| has discussed a whole sequence of higher order Lagrangians 

2" Wi^'i -^^U^Vi • • • -^^l^2rL-ll^2n\ V i/- K^'^' J 

This gives Ry/—g for n = 1 and reduces to ()29j) for n = 2. 

In general, /Z*^"^ vanishes identically for D <2n — \ and is a total derivative for D = 2n 
(the generalized Gauss-Bonnet theorem j3])- For D > 2n + 1, leads to second order 
equations of motion. The corresponding variations were also found by Lovelock 

5(1 rf^x) = I Gt'JSg'^'V^d^x, (31) 

/^W _ (2r2 + 1)! p V2n-lU2n ( On\ 

'^Xa ~ 2"+! y>^W-'^>^l>^2 ■ ■ ■ ^y2n-iy2n] ' K"^^ ) 

For n = 1 this coincides with the Einstein tensor. 

Explicit expressions for and G^^J {n < 4) with terms of equivalent tensor index 
structure collected as in can be found in j29i. We won't need them, since in practice 
for n > 2 it is much easier to work with (j3(J|) and (j32j) . 

Sometimes in the literature the theory (|29|) in D = 5 is attributed to Cornelius Lanczos, 
citing his papers 3^. In fact, however, Lanczos never went beyond D = 4.^ The name 
'Gauss-Bonnet' often associated with gravity theories based on Lagrangians (jHUj) is also an 
example of unfortunate terminology, since in D > 2n + 1 these 'dimensionally continued' 
Gauss-Bonnet densities are no longer associated with topological invariants (which is pre- 
cisely why they become interesting from the dynamical point of view). In this paper we will 
use the term Einstein-Lovelock gravity for the theories based on (j29|) . (j3(jp . 

It is often incorrectly stated that C^'^^ is the only combination of curvature squared correc- 
tions to the Einstein gravity consistent with string theory. This claim is based on a (correct) 



^ In 1932 he did not even include {R^i^xaY in Lagrangian; in 1938 he proved that 1)29(1 is a topological 
invariant density in _D = 4. I would like to thank N. Deruelle, J. Madore and J. Zanelli for the interesting 
discussions of this historical matter. 
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result of Zwiebach jSjJ, who showed that ()28p is the only curvature squared gravity theory 
in which the quadratic term in the expansion of the action around flat space is the same as 
in the pure Einstein gravity. In particular, the graviton is the only particle in the perturba- 
tive spectrum, and unphysical ghost poles usually associated with curvature squared terms 
(see, e.g., js^) are absent. However, one should remember that only on-shell effective action 
can be defined in string theory. This is especially clear when this action is computed from 
the S-matrix, but is also true for the sigma model /3-function method (see, e.g., |33]). Thus 
there always remains field-redefinition freedom, which can be used for example to change the 
coefficients in front of or completely remove {Rp.vY ^'^d -^^ terms in ()29j) . All these actions 
would be equally consistent with the string S-matrix and, thus, with string theory 3J]. In a 



sense, string theory does not have much to say about off-shell dynamics of quantum gravity. 



I thus prefer to think of higher curvature theories like (jzHj) not as fundamental microscopic 
theories, but as effective field theories of gravity (see, e.g., recent discussion in js^), with 
the hope that they may capture some aspects of gravitational field dynamics in presence of 
regions of high curvature. 



V. METRIC ANSATZ AND EQUATIONS OF MOTION 

We will now study the effect of the Lovelock term on the collision spacetime at 
M,f > 0. For simplicity, we will analyze the head-on collision (6 = 0). In principle, the 
method could also be used to study the nonzero impact parameter case, most certainly with 
similar conclusions. Remember that we always assume D > 5. 

Our collisions are transplanckian: E ^ 1. It is this situation that can lead to formation 
of a large classical black hole according to the common lore 0, S Q • For ~ 1 we would 
be speaking about Planck-size black holes, for which quantum gravity effects are significant 
without doubt. 

Our goal is to find the metric to the future of the collision plane up to u,v < 1, at 
transverse radii r ~ tq. Since r ^ 1, we can neglect transverse derivatives of the metric: 
at M,f ^ 1 they are suppressed by 1/r compared to the longitudinal ones [see (fT^ . (US))]. 
Neglecting transverse derivatives means that we approximate the metric near a given point 
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V in the collision plane located at r ~ ro by a plane wave collision metric (see j36[) 



D~2 



ds^ = e^("'")c/M dv- [1 + A{u, v)] {dy^f - [l + B{u, v)] ^(c/y')'. 



(33) 



i=2 



Here are Cartesian coordinates in the transverse plane near V, with the axis pointing 
in the direction of the collision point, and y"^, . . .|/^~^ in the D — 3 orthogonal directions 
(see Fig. El). 




FIG. 3: The transverse collision plane. Shown are the trajectories of the colliding particles and 
the circular intersection with the CTS ((T^ . 



The initial conditions to be imposed on L, A, B outside the future wedge F{u, f > 0) are 
obtained by matching (j33|l with (fT6|) . (fT5|l . We find {u,v ^ F) 



L{u, v) = 0, 

A{u, v) =e [ue{u) + v6{v)] , 



(34) 



B{u,v) 



D-3 



[ue{u) + ve{v)]^ 



where 



(35) 



^ [D - 3)E ^ i/(^_3) 
To derive the equations of motion, it is convenient to work with a more general metric^ 

D-2 

ds" = e^("'")rfM dv-Y, [^^'^ i% v)]\dy')\ (36) 



i=l 



^ In the rest of this section we suppress Einstein's convention about summing in repeated indices, and 
indicate all necessary summations explicitly. 
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The nonzero components of the Riemann tensor are {Au \ L™, . . . are partial derivatives) 

D _ 1 r 

R ■ ■ = M(*) - 4(*) L 1 A^^^ 

D . . _ A{i) A(.i) (37) 

p . . _ r /I (j) _ 4 («) r 1 4 (*) 

-ftjjjj — z e /i^ -t- /l^ /!„ J /i /I . 

We will now make an a priori assumption that 

|VL|, |VA(*)| < 1 {u,v< 1), (38) 

which also implies 

L^O, A^^^l {u,v<l). (39) 

Validity of this small- gradient approximation will have to be checked after a solution is found. 
For now we see that it is compatible with the initial conditions fl34p . 

The approximate Riemann tensor is found by applying (jHHj) and (IH^ in (jHIj). We have: 

R ~ — r 

R • . ~ /!(*) 

I? . . ~ /((*) 

B ■ ■ K. A^"^ 

all other components being 0. 

The vacuum equations of the Einstein-Lovelock gravity are 

+ KGf^ = 0. (41) 

The uv and ii components of this tensor equation will give D — 1 independent equations for 
D — 1 functions L and A'^'^^ . The nonzero components of Gao- and G^^l in the small-gradient 
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approximation are easy to find using (pUj). For Gao- we have: 

i 

'-J^uu — / ^ ^UUI 

(42) 

= -2L„, - 4 J] (1 < ^ < D - 2). 

("2") 

For G^-;^ we use the general expression (j32|) with n = 2. Because of the antisymmetrization 
in the r.h.s. of (jS2I), only terms for which cr, Ui,..., V2n are all different can contribute. 
Also, looking at (j4Uj) . we see that u or v must be present in each pair of indices (z/i, z/2), . . . , 
(^2n-i, V2n)- Usiug such reasouiug, it is easy to see that for n = 2 the only nonzero compo- 
nents are'^ 

gS? = 16 J] [A^A^ - {l<^<D-2). (43) 

We now go back to our original ansatz (|^. putting 

A(^) = 1 + A(u,?;), A(^) = 1 + 5(u,?;) (2 < i < D - 2). (44) 

Substituting (ji^ and in (jH)), we have three independent equations for L,A,B. The 
simplest is the ww-equation, to which only the Einstein tensor contributes: 

+ {D- 3)Bu, = 0. (45) 

This equation is easy to solve. Namely, it implies that the relation 

^/ X A(u,v) 

BM = -^±j^ (46) 

satisfied by the initial data (jH^ . will continue to hold for M,t> > 0. 

There remain two ii equations {i = 1,2). After B is excluded using ()46p. they can be 
brought to the form 

D — 3 

i^\_AuuAyy — (A^t,)^] + — 4y^"'' ~ ^' i'^'^) 
Luv = Auy. (48) 



Notice that for n > 3 the given conditions on indices are incompatible. This means that the higher order 
Lovelock tensors vanish in the small-gradient approximation. 
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After a function A{u,v) satisfying ()47|) is found, it will be easy to find L{u,v) from 
Taking the initial conditions into account, we will have 



L{u, v) = A{u, v)-e [ue{u) + ve{v)] . (49) 

Thus we reduced the problem to solving a single nonlinear second order PDE (|47p . In 
the pure Einstein theory {k = 0) the nonlinearity disappears, and we see that the solution 
for M, f > is given by the same formulas ()34|) as the initial data. For k ~ 1 this is no longer 
a valid solution, since it has 

Au^A^^ - Al^ = e^5{u)5{v), (50) 

which becomes large at u = v = (namely, ~ e'^E^ ^ 1, when the finite shock width 
w ^ 1/E and the resulting smearing of the 5-functions are taken into account). 

PDE fl47|) belongs to the well-known family of hyperbolic Monge-Ampere equations. It 
may be noticed (see 0]) that it can be derived from the action 

^MA = j ciurfw|| {AlA^^ - 2A^AuA 

= j dtdz[^e''''' e''' daAd,Ada>yA + -^^^^^,^-'daAd,A], (51) 

which we should consider as a 2-dimensional truncation of the full Einstein-Lovelock action 
relevant for our problem. The first term in this action uses only the totally antisymmetric 
tensor to form contractions. For this reason the corresponding nonlinear term in Eq. ()47p is 
invariant with respect to Euclidean rotations as well as Lorentz transformations. The linear 
wave equation term in (j47p is of course only Lorentz-invariant. 



VI. SOLVING THE MONGE-AMPERE EQUATION 



A. Smooth solutions 



Hyperbolic Monge-Ampere equations have exceptional character among the nonlinear 
second order hyperbolic PDEs: they can be reduced to a system of five first order charac- 
teristic ODEs, while in general (in two dimensions) one needs eight lad p. 495]. However, 
in case of Eq. (jTrf) . because it has constant coefficients, this general theory is superseded by 
an even simpler method based on the Legendre transform. 
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Let us rewrite (jTTj) as 



2 D — 3 

^uu-A-vv — ^™ + ^-^uv = 0, A = , -T— , (52) 

10(1^ — A)k 



substitute 



A = \{A + uv), (53) 
AuuAvv ~ — Auv = 0, (54) 

and apply the Legendre transform 0, p. 32]. This results in a linear Poisson equation, 
which is readily solved. In this way Ignatov and Poponin 0] obtained an exact solution of 
the original Monge- Ampere equation ()52|) in an implicit form, depending on two arbitrary 
functions -F(0 and G{ri): 

A{u, v) = X [F{0 + G{v) + F'iO G'iv)] , (55) 
where ^ and r] must be found from 

(56) 

v = ^ + G'iv)- 

For vanishing at infinity small functions F, G this solution describes interaction of two 
colliding pulses of the A field. The pulses emerge from the interaction region with their 
shape unchanged, resuming the motion along the pre-coUision world lines. 

In general, solution ()55|1 is applicable if the transformation (?7, ^ {u,v) given by ()5fi|l 
is one-to-one. This means that its Jacobian should not vanish: 

|^ = l-F"(e)G"(r,)^0. (57) 

Unfortunately, this condition is violated in our problem. Indeed, to satisfy the initial con- 
ditions we would have to take 

F{0 = G{0= 6^6(0. (58) 

Because of finite shock width w ^ 1/E these initial data have to be smeared out on this 
scale. Still this results in F"{C,) near ^ = of the order ~ eE ^ 1. So the Jacobian 
()57|) . being equal to unity away from the origin, changes sign and becomes negative near 
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^ = t] = 0. The transformation ()56p is thus not one-to-one: there is a region of the u, v 
plane which is covered 3 times, and where the inverse transformation is muhiple-valued. 

This analysis shows that for initial data (|34j) . even after smearing them out, the exact 
solution (j55|) in any case cannot be used verbatim. A region of spacetime appears where 
this solution is triple-valued, and we have to somehow choose between the three branches. 
Actually, it turns out that no satisfactory choice is possible. More precisely, any choice 
introduces discontinuities either in the function A itself (which is totally unacceptable) or 
in its first derivatives (which brings back the problems which we temporarily resolved by 
using smoothed out initial data). All of them seem to violate Eq. ()52|1 . 

To get some insight about the source of these difficulties, I also studied the problem 
numerically. The conclusion of this study (which I am not going to report here in detail) is 
that the solution develops singularities characterized by discontinuous first partial deriva- 
tives. Evolution beyond these singularities apparently cannot be described by formula (|^. 
A different strategy is required, which I will now proceed to describe. 



B. Generalized solutions 



To find a solution, we will have to extend the class of admissible functions A, allowing 
continuous functions with discontinuous first derivatives. [Notice that the initial data ()34|) 
corresponding to infinitely thin shocks belong to precisely this class.] Such generalized solu- 
tions, called weak in the theory of PDEs j3, P- 418, 486], are usually defined via integration 
by parts. In our case weak solutions can be defined, because the Monge-Ampere operator 
can be written as a divergence 

Auu^vv ~ ~ dv{AyAuu) ~ du{AyAuv)- (59) 
Thus Eq. ()52p is equivalent to the conservation law 

duJy + dyJu = 0, (60) 
where the current has components 

Jv — ~AyAuy + ^^r>) (^l) 

Jy = AyAyy + ^A^. (62) 
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For smooth functions A the differential form ()60|) of the conservation law is equivalent to 
the integral form: 

^e'^Varfxfe = (63) 

for any closed curve C. However, the latter condition does not involve products of second 
derivatives of A and can be unambiguously checked even when the r.h.s. of (j5!2j) cannot be 
defined. I thus take as the definition of weak solutions of Eq. fl52|) . 

In practice, it is always sufficient to check ()63p only for infinitesimally small contours 
C, and this only near special, most dangerous, points. In the rest of the u, v plane the 
differential equation (|52j) can still be used. Provided that the first partial derivatives of A 
are bounded, the nonsingular terms proportional to A drop out of (j63|) in the limit of small 
contours. The remaining rule for catching possible 6{u) )-type singularities as in (|5Up 
takes the form 

J {AuuAy^-Al^)dudv= j AydAu. (64) 

[This formula has in fact general validity and follows directly from (jSHl)-] If there is a 6- 
function at the origin present in the integrand on the l.h.s., we will be able to detect it by 
studying the limit of the r.h.s. as e — 0. The latter procedure does not involve products of 
distributions and is unambiguous. 

For further discussion it is convenient to simplify ()52p by making the substitution 

A = A{u,v) + ^uv, (65) 
AuuA.. - Al, + AV4 = 0. (66) 

This PDE is hyperbolic due to positivity of 

We will no longer attempt to use smeared out initial data, and instead will try to find a 
solution directly in the limit of infinitely thin shocks. Since in this limit u = v = is likely 
to be the most singular point, we change to polar coordinates 

u = psiinp, V = pcosip. (67) 

PDE becomes 

^ [Ap(^w + pA)- ( - ^^/pY] + ^74 = 0. (68) 
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The crucial fact (easy to guess if one remembers the well-known connection between the 
homogeneous Monge-Ampere equation and developable surfaces) is that for functions A of 
the form pf{(f) the first most singular term in (j68|) vanishes. To satisfy the full equation, 
we have to include in A terms of order p^. Thus we are led to the ansatz*^: 



Compatibility with the initial conditions (j34j) requires: 



(69) 



e sin yj, G [vr/2, vr], 
0, G [7r,37r/2], 

ecos(/9, if G [37r/2,27r] 



(70) 



A 



sin2(/9, if G [vr/2, 27r]. 



For If G [0, 7r/2] the functions /, g have to be found so that PDE is satisfied (see Fig.Hj). 

/Ml 



? 



7r/2 



37r/2 



27r 



<y3 



FIG. 4: The initial data for f{'^). In the future wedge F (see Fig. ^ the function has to be 
determined (?). The dashed line denotes the Einsteinian solution (jTH]). 



Checking PDE (pj)) at m = w = calls for the use of identity An easy calculation 

shows that 

(71) 



It is possible to consider a more general ansatz, adding terms p"ft,„ ((/?), n > 3. The functions hn would 
then be restricted by an infinite sequence of nonlinear ODEs, obtained similarly to 1)74(1 and ((75|l below. 
These higher order terms can be neglected compared to ((69|l in the limit of small p, and so we put them 
to zero for simplicity. In principle, however, the possibility of them being nonzero further aggravates the 
non-uniqueness issue to be discussed in Sec. VI. C below. 
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Taking (fTUj) into account, we get a condition for the absence of S{u)S{v) in 

r^\f"-f)dv = 0. (72) 
Jo 

As an example when 6{u)6{v) does appear, we can take the Einsteinian solution A = 
e[u9{u) + v9{v)~\, which would correspond to 

f{(p) = e {simp + cos (p), (pe[0,n/2]. (73) 

And indeed we readily see that (f?^ is violated. 

At p > the solution can be understood in the usual sense, and the use of is 
unnecessary. We simply plug the ansatz (fH^ll in Eq. (jHH|l . Equating the coefficients before 
p^^ and to zero, we get two ODEs: 

(/" + f)9 = 0, (74) 
2gg"-g" + Ag' + XyA = 0, (75) 

which have to be solved with the following from (|7Up boundary conditions 

/(O) = /(7r/2) = e, (76) 
g{0) = g{n/2) = 0. (77) 

At the first glance we have a problem here, since ODE ()74|) implies that either g = ot 

f" + f = 0. (78) 

The first opportunity seems to be excluded by (fTSjl . while the only solution of (fTHj) with the 
boundary conditions (ffHjl is (ffHjl . which does not satisfy the no-delta condition (ff^ . 

The way out of this impasse is to notice that although g is indeed not allowed to be 
zero on an interval, it can still vanish at isolated points inside (0, vr/2). At these points the 
derivative f'{f) can have jump discontinuities. If on the rest of the interval (fTSj) holds, ODE 
dZH) will still be satisfied. 

The general solution of ODE (jTSj) depends on two constants C, (f and has the form 

g{ip) = ±^{VC^cos[2{if - ^)] - C}. (79) 
At points ip^: where g{ip^:) = we have, in agreement with (fTSj) . 

g'{ip,) = ±A/2. (80) 
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The constants C, (p and the overall sign in ()79|) can be chosen independently on the two 
sides of provided that both choices are compatible with (7(^2*) = 0. In particular ^''(v?*) 
is allowed to flip sign at <^*. Being multiplied by vanishing the arising 5- function in g" 
does not violate (f7H|) . 

Taking all these observations into account, we can write down the following class of 
solutions to the boundary value problem ()74|) - (f77j) : 

• Functions /, (7 are continuous on [0,7r/2]. 

• The / satisfies f|75|) everywhere on [0,7r/2] except at some chosen points 

< V9i < . . . < v^AT < 7r/2, (81) 
where will have jump discontinuities. 

• The g vanishes at y^i, . . . , (/^at, and is given by (f?^ on each subinterval v^a+i] 
(v^o = 0, v^AT+i = 7r/2 are the end points) with 



(82) 
(83) 



C = COt{(pa+l - (fa), 

and an arbitrary overall sign ±, which can be chosen independently on each subinterval. 

/Mr 





-"^^^K = 1 + a/2 




K = 1 - a/2 


Fi 


F-2 



7r/2 



FIG. 5: Solution (jSlJ). The future wedge F is split into two equal sectors Fi^2 (see Fig. P). The 
leading behavior of A{u, v) in these sectors is given hy u + Kv and v + Ku. 

Within this class it becomes possible to satisfy the no-delta condition (f7^. The simplest 
and most symmetric example can be constructed using one intermediate point ipi = it /A 
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(see Fig. EI). From continuity and ()78p we must have 



cos ip + K simp, (pE [0,11/4], 



(84) 



sinip + K cosif, G [7r/4, 7r/2]. 
The constant K is determined from (f?^ . We find that two values are allowed: 



K = 1±V2. 



(85) 



For gl^p) we have: 



A I v^cos(2(p-7r/4)-l, y;G[0,7r/4], 
^ ' y2cos(2(^ - 37r/4) - 1, (p e [tt/A,tt/2]. 



(86) 



Near the point u = v = the leading pf{(p) behavior of these solutions (as of all the solutions 
in the class described above) is piecewise linear in terms of the u, v coordinates, with the 
p'^g{(p) term being a small correction (see Fig. IHI). 



Einstein 




K = ^ + ^/2: 



FIG. 6: The leading behavior of A{u, v) near u = v = Q for the pure Einstein case (|73|) and for the 
Einstein-Lovelock solution 1)84(1 . The time flow is from the left to the right. 
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A more general example, which we will need in a discussion below, utilizes two symmetric 
intermediate points y^i = 0, (/92 = 7r/2 — 0, where G [0, 7r/4]. The function / is given by 



f{if) = ex < 



cosy? + i^sinyj, 
cos d + K sin 6 



if e [0,0], 

cos ip + simp) , G [0, 7r/2 — 0], 



cos + sm ( 

simp + K cos ip, V9 G [7r/2 — 0, 7r/2], 



(87) 



where the K is fixed by (f?^ to be 



K =1± v^l + coti 



(88) 

7r/4 



The allowed functions g can be written down according to the rules given above. For < 
this example reduces to the previous one. 

Taking even more general configurations of intermediate points and imposing the no- 
delta condition ()72j). we can get infinitely many solutions of PDE ()66|) having the form of 
the ansatz and satisfying the initial conditions (f7n|) . Solutions of the original Monge- 
Ampere equation (fS^ are then given by Eq. 



C. Discussion of non- uniqueness 

Let us repeat the logic of the preceding discussion. In Sec. VI. A we convinced ourselves 
that solutions of the Monge-Ampere equation (j52|) with the smeared out initial conditions 
(jH^ run into singularities, characterized by discontinuous first partial derivatives. The very 
meaning of a solution in presence of such discontinuities needs to be reconsidered. For 
this reason in Sec. VLB we introduced the definition of weak solutions, using a standard 
procedure of integration by parts. ^ Working in the limit of infinitely thin shocks, we found a 
whole class of weak solutions with the required initial conditions. It is of course satisfying to 
find that solutions according to the new definition exist, but how should we interpret their 
non- uniqueness? 

In my opinion, there are two alternative possibilities. First, it is possible that this non- 
uniqueness can be resolved by performing a detailed study of PDE ()52j) with smeared out 

^ Notice that for smooth functions the concepts of usual and weak solutions are equivalent. It is only in 
presence of discontinuous first derivatives that the usual notion is not applicable, and we have to use the 
new definition. 
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initial data at small but finite shock width w ^ 1/E. As we have just discussed, the 
singularities will still appear. But perhaps by studying their approach one may discover a 
preferred way to continue the solution past them. This would require techniques beyond 
what we use in this paper. However, we would like to point out that even if a 'true' unique 
solution is found in this way, it will still have to satisfy our weak solution criteria from 
Sec. VLB. Because of this it seems likely (although cannot be guaranteed) that this 'true' 
solution will be close to one of the weak solutions we found in Sec. VLB. 

The second, more intriguing, possibility is that the encountered non-uniqueness is funda- 
mental. Then its interpretation can be most naturally given in quantum theory. It would 
simply mean that beyond the shock collision the gravitational field wavefunction ceases to 
be concentrated on a single classical configuration. Instead, it spreads out over all allowed 
solutions of the equations of motion. In this case we would have to sum over all solutions, 
using exp(iS'), where S is the on-shell value of the action, as the weight. 

It should be noted that the problem of non-uniqueness of weak solutions is well known 
in the theory of PDEs. It is present for instance in the theory of first order quasilinear 
hyperbolic systems, where the weak solutions in question are shock waves. There the non- 
uniqueness is resolved by imposing the physically motivated 'entropy condition' |4o| . 

One can try to do something similar in spirit for our PDE (jH21)- For example, one may try 
to identify higher order currents that are conserved for usual smooth solutions, and impose 
their conservation as an additional constraint on weak solutions. In fact, the starting point 
of our discussion of weak solutions was to rewrite the Monge- Ampere equation ()52j) itself as 
the current conservation condition Going a step further, let us introduce a symmetric 
2-tensor Tab with components 

Tuu = (Vv4) Auu, 

= {VA)'A,,, (89) 



Al), (90) 
-Al). (91) 

Actually, the appearance of the Monge-Ampere operator in the r.h.s. is not surprising: 
Tab is nothing but the energy-momentum (EM) tensor of the homogeneous Monge-Ampere 
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We have 



Tuv — T^u — ~(^A)'^Ajj 

^v^uu ~t~ ^u^vu '^Ay(^Aij^'aAijy — 



equation, Eq. ()52|) with A = 0. This equation corresponds to the first term in the action 
fj51|) . To define the symmetric EM tensor, we covariantize this term: 

j d^xe'"''e^^'daAdi,Ada'yA^ j d?Xy/^ E"^' E^^' daAd^AV a'dyA (92) 

and differentiate with respect to the auxihary 2-dimensional metric Qah- The result of this 
computation is (jHI?|) (up to a factor of 3/2). 

The EM tensor of the full equation (|52|) differs from (|89p by a trivial harmless piece 
associated with the quadratic term {dAY in (jHT|l . In imposing the EM conservation at the 
singular point u = v = Q that piece is irrelevant. Transforming the conservation laws (jHOl), 
()91|) to their integral form, we conclude that the EM tensor will be conserved if 



{VAydAu, f {VAydA^^O (e^O). (93) 

Applied to the ansatz (jU^ . this gives two conditions on f{ip), conveniently written down as 
one complex condition: 

^d^e^nf' + n-if + n=0. (94) 

When this condition is imposed, the class of weak solutions from Sec. VLB is reduced, but 
not eliminated. The simplest solution ()84|1 does not satisfy ()94j) and would not be allowed. 
However, already among the one-parameter family (j87p there is a solution [corresponding to 
the numerical value (p ~ 0.322 and the minus sign in (j88p ] which passes the new criterion. 
Introducing more intermediate points, we will still have an infinitude of solutions satisfying 
both the no-delta condition (f7^ and the EM conservation ()94|) . 

It seems likely that a whole hierarchy of currents conserved on the smooth solutions of 
PDE fjH^ can be found. In fact, it is possible ||41|] to map the Monge-Ampere equation 
to the Born-Infeld equation 

(1 + ^D^u - 2^t^,^t. - (1 - $?)$xx = 0, (95) 
or to the first order system describing dynamics of the Chaplygin gas 

Pt = PPx + Q ^Qx, 

(96) 

Qt = {PQ)x- 



Both (j95j) and (j96|) are known to possess infinitely many conserved quantities j42]. Whether 
translating them back to the Monge-Ampere variables and imposing on the weak solutions 
is physically motivated, or can help resolve the non-uniqueness, is an open question. 
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This is as far as our explorations have taken us in the attempts to resolve the issue 
of non-uniqueness. Unfortunately, the problem still remains. In particular, at present we 
are unable to decide which one of the two possibilities outlined at the beginning of this 
subsection is realized. However, both of them allow for comparison with the pure Einstein 
case, to which we now proceed. 

VII. APPROXIMATE POST-COLLISION METRICS COMPARISON 

In Sec. V, when deriving the equations of motion, we introduced a small-gradient as- 
sumption (|38|). which for our particular ansatz (|33|) takes the form 

|VA|, |V5|, |VL| < 1 iu,v<l). (97) 

This was compatible with the initial data, but could not be taken for granted for the full 
solution. So before we use the solutions found in Sec. VLB, we have to check ()97|) . 

The functions B and L are trivially related to A by ()46p and ()49p. Thus it suffices to 
check fl97|) only for A{u,v), which is given by ()65|) and (jSHI)- To estimate the gradient of 
(jnni), we compute 

V4p/(^)] = fiv) sin^ + f\4>) cos<^, (98) 
V« [p^g{v)\ = p [25'(v?) sin ip + g'{(p) cos (p] (99) 

Analogous formulas with sirup and cosip interchanged are valid for V^,. Since f{<p) will 
typically be of order e <^ 1 [as in (jS^ and (jS7j) ]. the pf{(p) term passes the small-gradient 
check. Further, the size of giip) is set by the parameter A, which will be ~ 1 for k ~ 1. It 
follows that the p'^gi'p) term, as well as the term Xuv/2 in ()65|) . will also pass the check, 
provided that we restrict u,v to a. somewhat smaller region 

M,W^Pmax<l, (100) 

where the precise value of pmax has to be chosen depending on the required accuracy. This 
value is insensitive to E, which determines the shock width w ~ 1/E. This means that 
the found solutions do take us out of the shock interaction region before the small-gradient 
approximation breaks down. 

The plane wave collision metrics ()33p corresponding to the solutions of Sec. VLB approxi- 
mate the spacetime near the chosen point V in the transverse plane. They implicitly depend 
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on V through the parameter e. To assure the consistency of neglecting transverse deriva- 
tives, we have to assume that all other free parameters, such as the choice of intermediate 
points (fai depend weakly on or are independent of V. Then we can patch the solutions for 
different P's together. The resulting approximate metric has the form 

ds" = e^("'"''^)du dv- [I + A{u, V, r)] - [l + B{u, v, r)] (101) 

where the dependence of L,A,B on the transverse radius is now made explicit. These are 
the metrics described in the Introduction. The final range of their validity is 

r~ro, M,W^Pmax- (102) 

In the pure Einstein case, the approximate post- collision metric at r ~ tq, f ^ 1 would 
have the same form ()101|) . except that the function A has to be found from the linear wave 
equation A^v = (see Sec. V). Thus we have 

Ae{u, V, r) = e{r) ■ [ue{u) + ve{v)], (103) 
A-F(u,v,r) 

BE{u,v,r) = (104) 

LE{u,v,r) = 0. (105) 

Comparing the pure Einstein and the Einstein-Lovelock spacetimes, we see that they are 
quite different. In the Einstein case the Riemann tensor is concentrated on the null planes, 
while in the Einstein-Lovelock case it has additional 5-function singularities of comparable 
strength on the planes u/v = tampa corresponding to the intermediate points (jHH)- This 
difference persists as far as we can trace the solutions, that is up to M,f ~ Pmax- 

It is instructive to see what happens with the dynamics of CTS in such a spacetime. 
In Einstein gravity, once a CTS appeared, it can be evolved into the future along its null 
normals. The CTS property is preserved by such evolution because of the Raychaudhuri 
equation (j^H). This fact crucially depends on the inequality 

Rt^J'^i" > (106) 

(for any lightlike £^), which follows from Einstein's equations when the matter sources are 
either absent or satisfy appropriate positivity conditions jlfi| . 
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The Ricci tensor of spacetime ()10ip is best written down in the same local coordinates 
u,v,y'' that we used in Eq. (jSSl)- Neglecting transverse derivatives, and using the small- 
gradient approximation, we have 

F? ~ — r — — A 

Rylyi ^ AA^^, (107) 
4A 

i?,V~45™ = --^-^ i2<t<D-2) 

[where we used fl46|) and to express B and L via A]. We see that in the Einstein- 
Lovelock case R^i, 7^ in the future wedge. Most notably, R^i, will have singularities on 
the planes u/v = taiupa, with the leading behavior oc p5{ip). The magnitudes and signs of 
these singularities depend on the particular choice of solution. 

To see what happens with condition (jlOfjj) . let us pick a light-like vector ={i'^,i'",t). 
Since the metric components do not differ much from their Minkowskian values (only their 
second derivatives do), we have 

rr - (f )2 ^ 0. (108) 

Using p07|) and ()108|) . it is easy to compute 

r D - 1 ^"^ 1 

R,J^r = 2Au, [{ey - ■ (109) 

i=2 

The signs of both factors in this expression are undetermined, and so Ricci tensor ()107|1 will 
generally violate condition pOfjj) . 

We thus see that in the considered example the higher curvature terms significantly 
modify the post-collision spacetime. It appears that the pure Einstein gravity and the CTS 
argument based on it cannot be trusted in the future of the collision plane. 



VIII. FINAL REMARKS 



In this paper I questioned the use of classical Einstein gravity to analyze black hole 
production processes in transplanckian particle collisions. My basic argument was very 
simple: I pointed out that curvature, as measured by the curvature invariant {R^uXaY, 
becomes large on the transverse plane u = v = in the classical collision spacetime, in the 
region relevant for the horizon formation. 
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In principle, I could have stopped right here, since this fact alone already implies that 
quantum corrections cannot be ignored. However, I decided to go a bit further and give a 
concrete example of how taking quantum corrections into account may dramatically change 
the picture. 

Needless to say, the precise form of quantum corrections to classical general relativity 
remains so far unknown. However, there is a common opinion (supported by string theory) 
that at least to some extent these corrections may be represented by adding to the effective 
action of gravity a sequence of higher curvature terms. 

In this paper we considered a particular example of such 'higher curvature gravity' — the 
so-called Einstein-Lovelock theory. We re-analyzed the collision in this theory and found 
that the spacetime does deviate significantly from the pure Einstein case. Thus, the initial 
conclusions are supported. 

Of course, once it is found that a certain combination of higher curvature terms strongly 
modifies the Einstein gravity result. Pandora's box is open, and there is no reason to expect 
that all other higher curvature terms can be excluded from the analysis. Moreover, since 
the width of the colliding shock waves is ~ E-^ < 1, it seems very likely that the higher 
the power of curvature, the more important this term will become. (It was not so for the 
higher Lovelock Lagrangians solely because of antisymmetrization.) 

For this reason, the given analysis of the Einstein-Lovelock gravity case should not per- 
haps be assigned importance beyond that of an example supporting the main argument. 
(Nevertheless, it seems to be the first study of a dynamical process in this theory; as far 
as we know, all previous applications involved static solutions.) The real conclusion of this 
analysis and the whole paper is that, at the present state of knowledge about quantum 
gravity, no well-founded claims about transplanckian collisions at small impact parameters 
can be made. Even if TeV-scale gravity is realized in nature, the true story of such collisions 
will be much more complicated than a single large black hole production with a subsequent 
Hawking evaporation. 

Perhaps some readers will find the following additional explanation useful. It is sometimes stated that 
quantum gravity corrections in transplanckian collisions should be small because the Schwarzschild radius 
-Rschw(-E) » 1 for £^ 3> 1. This would indeed be true were Rschw{E) the only relevant length scale in the 
problem. But my analysis shows that it is not: the other relevant scale is set by the curvature radius near 
u = t; = Oatr~ Rschw{E), which according to (j?7j) is actually much smaller than the Planck length. 
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